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We use the Snla Gold dataset to compare LCDM with 10 representative parametrizations of the 
recent Hubble expansion history H(z). For the comparison we use two statistical tests; the usual 
Xmin an d a statistic we call the p-test which depends on both the value of Xrnin and the number n of 
the parametrization parameters. The p-test measures the confidence level to which the parameter 
values corresponding to LCDM are excluded from the viewpoint of the parametrization tested. For 
example, for a linear equation of state parametrization w(z) = wo + wiz the LCDM parameter 
values (wo = —1, wi = 0) are excluded at 75% confidence level. We use a flat prior and Dom = 0.3. 
All parametrizations tested are consistent with the Gold dataset at their best fit. According to 
both statistical tests, the worst fits among the 10 parametrizations, correspond to the Chaplygin 
gas, the brane world and the Cardassian parametrizations. The best fit is achieved by oscillating 
parametrizations which can exclude the parameter values corresponding to LCDM at 85% confidence 
level. Even though this level of significance does not provide a statistically significant exclusion of 
LCDM (it is less than 2a) and does not by itself constitute conclusive evidence for oscillations in the 
cosmological expansion, when combined with similar independent recent evidence for oscillations 
coming from the CMB and matter power spectra it becomes an issue worth of further investigation. 

PACS numbers: 


I. INTRODUCTION 

Converging observational evidence that appeared dur¬ 
ing the past decade has indicated that we live in a 
spatially flat universe with low matter density that 
is currently undergoing accelerated cosmic expansion 
[1, 2|@, 4s B 6,Q. This accelerating expansion has been 
attributed to a dark energy H] component with nega¬ 
tive pressure which can induce repulsive gravity and thus 
cause accelerated expansion. 

The simplest and most obvious candidate for this dark 
energy is the cosmological constant A (9j with equation of 
state w = p/p = — 1 . Such a model predicts an expansion 
history of the universe which is described by an expansion 
rate H(z) as a function of the redshift z given by 

H"(z; Ulom) = = •®o[^0m(l + z ) 3 + (1 — f^Om)] 

( 1 . 1 ) 

where flatness has been imposed and Dom = po/pc is 
the single free parameter of this simplest data consistent 
parametrization (LCDM). 

The most sensitive observational probe for testing this 
type of parametrizations comes from distant standard 
candles like type la supernovae (Snla). These make it 
possible to start seeing the varied effects of the universe’s 
expansion history. The most updated and reliable com¬ 
pilation of Snla is the Gold dataset recently relased by 
Riess et al. 0. The authors have compiled a catalog 
containing 157 Snla with z in the range (0.01,1.75) and 
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visual absorption Ay < 0.5. The distance modulus of 
each object has been evaluated by using a set of cali¬ 
brated methods so that the sample is homogenous in the 
sense that all the SNela have been re-analyzed using the 
same technique. Thus the resulting Hubble diagram is 
indeed reliable and accurate. Even though LCDM pro¬ 
vides the simplest parametrization consistent with the 
Gold dataset it has two disadvantages which motivate 
the search for other models: 

• It requires extreme fine tuning of the value of the 
cosmological constant A (coincidence problem). 

• It does not provide the best possible fit to the Gold 
dataset. 

In an effort to address these disadvantages three ap¬ 
proaches have been followed 

• Assume the existence of a homogeneous time de¬ 
pendent scalar field whose dynamics is determined 
by a specially designed potential so that its energy 
comes to dominate at present and its negative pres¬ 
sure plays the role of dark energy. Such fields could 
either have w > —1 if the sign of the kinetic term is 
positive (quintessence £ifj) or w < — 1 if the sign of 
the kinetic term is negative (phantom fields llll|). 

• Consider extensions of general relativity mill in 
mill (motivated eg by extra dimensions) which 
can lead to an accelerating expansion without any 
modification of the energy momentum tensor. 

• Examine arbitrary parametrizations mm of the 
expansion history H(z) and focus on maximizing 
the quality of fit to the Gold dataset. 
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The problem with the first two approaches which exhibit 
some physical motivation is that despite their increased 
number of parameters they very rarely exceed the qual¬ 
ity of fit of LCDM. This is not the case for the third ap¬ 
proach. Parametrizations which are not constrained to 
emerge from currently known physical theories provide 
significantly better fits to the Gold dataset than LCDM. 
For example it is very hard 0 to construct a physically 
motivated model predicting an evolution of the dark en¬ 
ergy equation of state w(z) that crosses the line w = — 1 
(phantom divide line). Yet the fit to the Gold dataset im¬ 
proves significantly if such crossing is allowed. The origin 
of this problem could either be statistical or physical. In 
the former case we are dealing with a statistical fluctu¬ 
ation of the data and upcoming more detailed datasets 
will become more consistent with some of the current 
theories (either minimally coupled scalar fields or modi¬ 
fied gravity). In the later case however, if eg the w = — 1 
crossing is verified by more detailed data, new theories 
will have to be constructed. Theoretical attempts in that 
direction have already started developing using either less 
physically motivated approaches involving combinations 
of scalar fields with positive and negative kinetic terms 
(2(1 21] or reconstruction of scalar tensor theories of grav¬ 
ity [IHlH 

In an effort to forecast and support this development 
it is important to extract the maximal information from 
the Gold dataset. In particular we wish to address 
the following question: ‘Given the number of parame¬ 
ters what is an arbitrary (flat and low matter density) 
H(z) parametrization that provides the best fit to the 
Gold dataset?’ Previous studies attempting to address 
this type of questions 00 have either been based 
0] on earlier less reliable datasets [2iJ [zij or have not 
been as extensive in the search of parametrization space 
0 Here we show that a careful and extensive search of 
parametrization space can reveal interesting new features 
of the expansion history hidden in the Gold dataset. 

We continue and extend previous work by two of us 
0 ] and evaluate the quality of fit to the Gold dataset 
for 10 representative parametrizations of H(z). They 
include physically motivated parametrizations like brane 
world 0. Chaplygin gas (first considered in [25), see 
also |26L l27j l and Cardassian cosmology, 0] as well as 
arbitrary parametrizations like a linear fit to the equation 
of state ( w(z ) = wq + w\z). 

The parametrizations compared have a different num¬ 
ber of parameters and therefore their comparison on the 
basis of how low a Xmin they achieve would be biased 
towards parametrizations with larger number of param¬ 
eters. The reduced Xmin (Xmin P er degree of freedom) 
is an improved statistic in that it has some weak depen¬ 
dence on the number of parameters. We propose however 
an alternative statistic, we call ‘p-test’ which avoids the 
above bias and is more appropriate for the comparison of 
parametrizations with different number of parameters. 
This statistic along with the general method used for 
evaluating the quality of fit to the Gold dataset of the 


parametrizations considered is described in the next sec¬ 
tion. In section III we present our results for the best fits 
and discuss their implications and common features. 


II. FITTING PARAMETRIZATIONS TO THE 
GOLD DATASET 


Given a parametrization H{z ; a ±,..., a n ) depending on 
n parameters we can obtain the corresponding Hubble 
free luminosity distance 


D tb (z;cn, ...,a n ) = (1 + z) [ dz' 

Jo 


Hn 


H(z'-,ai,...,a n ) 


( 2 . 1 ) 


Using the maximum likelihood technique 0] we can find 
the goodness of fit to the corresponding observed D°^ s {zi) 
(i = 1,..., 157) coming from the Snla of the Gold dataset 
jlj. The observational data of the gold dataset are pre¬ 
sented as the apparent magnitudes m(z) of the Snla with 
the corresponding redshifts z and lcr errors The 

apparent magnitude is connected to Dl(z) as 


m{z\ oq,..., a n ) = + 5log w {D L (z; oi,..., a n )) 

( 2 - 2 ) 

where M is the magnitude zero point offset and depends 
on the absolute magnitude M and on the present Hubble 
parameter Hq as 


c H -1 

M = M + 5log w (—M + 25 (2-3) 

Mpc 

The goodness of fit corresponding to any set of parame¬ 
ters ai,...,a n is determined by the probability distribu¬ 
tion of ai,..., a n i.e. 

P(M , Oi, ..., On) = (2.4) 


where 


157 

X 2 (M,ai,...,«„) = ^2 
2=1 


(m obs {zi) - m th (zi ; M, on,..., a n )) 2 
2 

a m° bs {zi) 


( 2-5) 

and A/" is a normalization factor. If prior information 
is known on some of the parameters ai,...,a n then we 
can either fix the known parameters using the prior in¬ 
formation or ‘marginalize’, i.e. average the probability 
distribution around the known value of the parame¬ 
ters with an appropriate ‘prior’ probability distribution. 
The parameters a \,..., a n that minimize the y 2 expression 
(1231) are the most probable parameter values (the ’best 
fit’) and the corresponding x 2 (ai,..., a n ) = Xmin gi yes an 
indication of the quality of fit for the given parametriza¬ 
tion: the smaller Xmin the better the parametrization. 

The minimization with respect to the parameter M can 
be made trivially by expanding |30) the x 2 of equation 
(1231) with respect to M as 


X 2 (ai, .., a n ) = A — 2MB + M 2 C 


( 2 . 6 ) 
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where 


^4 (m obs (zi) -m th (zi\M = 0,ai,..,a n )) 2 

A{a = 2 ^- 


P(ai,..,a„) = ^2 


i=l u m ohs {zi) 

(' m obs (zi ) - m th {zi\M = 0,ai,..,a„)) 


i= 1 
157 


m obs (z ■ 


c = 

i=l a m ob3 (zi) 

Equation lt2~6l) has a minimum for M = B/C at 


(2.7) 


X 2 (ai, a n ) = A(a i,a„) - (2.8) 

Thus instead of minimizing x 2 {M, ai, a n ) we can mini¬ 
mize % 2 (ai, a n ) which is independent of M. Obviously 
Xmin = Xmin- Alternatively we could have marginalized 
over the nuisance parameter M thus obtaining 1,^. Fill. ) 12 


X 2 (ai, ...,a n ) = ^ ^’ f \v,iCi'l~) 

(2.9) 

to be minimized with respect to Oi, a n . In our analysis 
we consider the y 2 (ai, ...,a n ) of equation (12.811 which is 
already minimized with respect to M. In what follows 
we omit the symbol tilde (7) for simplicity. 

Clearly the value of Xmin depends on the parametriza- 
tion used, and it can decrease arbitrarily close to 0 for 
large enough number of parameters n. Therefore, the 
value of Xmin by itself is not a particularly useful measure 
of the quality of a given parametrization. It is biased to¬ 
ward parametrizations with larger number of parameters 
and is therefore useful only in comparing parametriza¬ 
tions with the same number of parameters n. Notice 
however that the reduced Xmin ( Xmin P er degree of free¬ 
dom) has some weak dependence on the number of pa¬ 
rameters and can be an alternative when comparing the 
quality of fit of parametrizations with different number 
of parameters. 

How can we compare parametrizations with differ¬ 
ent number of parameters in a more efficient way? 
This can be achieved by the following procedure: Con¬ 
sider two parametrizations of the expansion history 
with different number of parameters a \,..., a n ) and 

H 2 (z-,b 1 ,...,b m ) and assume that for the parameter val¬ 
ues af,..., and ..., 6^ both parametrizations reduce 
to a common form eg LCDM i.e. 

P 2 (*;af,...,a£) = fff (*;&£, ...,&£) = (2.10) 

= ^Om) — Hq [flom(l + Z ) 3 + (1 — flom)] 

Let also Xirnin an d ximin be the minimum values of 
X 2 for each parametrization and y 2 the value of x 2 cor¬ 
responding to LCDM. Since Hl has a single parameter 
(flom) which is constrained by prior information (other 
observations) it is natural to expect that Ximin an d X 2 min 
will be smaller than y|. Define also Ay 2 = Xi — Ximin ■ 
A xl = xl - ximin > a Xil = xl - Ximin and A xl L = 
Xl — ximin • Now Ay 2 and Ay| are random variables 


that obey a x 2 probability distribution with n and m de¬ 
grees of freedom respectively. Therefore the probability 
that Ay? does not exceed the value A\f L is 


Pi<(&XiL,n) = 


1 


rAx?z,/2 


= 1 - 


2T(n/2) J 0 

T(n/2,A X IJ2) 


(f) 


l-i 


T(n/2) 


2 dx = 


( 2 . 11 ) 


and similarly for AyS;. In eq. (12.1111 T(n/2, Ayf L /2) is 
the incomplete T function. For n = 1 and Ay^ = 1 we 
find Pi<(l, 1) = 0.683 corresponding to lcr as expected 
from well known tables |2f)]. 

Therefore Pi<(Ay 2 L , n) expresses the probability that 
given the n parameter parametrization Hi(z\ a\, ..., a n ), 
the true values of parameters oi ,...,a n produce a Xi 
that is less than the value of y \ L corresponding to 
LCDM and therefore LCDM is not realized in na¬ 
ture. In other words it is the confidence level to 
which LCDM is excluded, from the viewpoint of the 
parametrization H\{z\ ai ,..., a n ). It is easy to see that 
Pi < (Axf i ,n) increases with Axf L but decreases with 
n i.e. a parametrization with many parameters has 
more difficulty to exclude LCDM unless it can provide 
a very low xLm ( lar g e A xIl)- Thus p i<( A xl L , n ) Pro¬ 
vides a quantitative measure of the quality of a given 
parametrization H\(z\ ai,..., a n ) which allows compari¬ 
son with a different parametrization H 2 (z] b i,..., b m ). We 
call this measure ‘p-test’ and it will be used in what fol¬ 
lows to compare the quality of the parametrizations con¬ 
sidered. 

Another way of comparison of the cosmological models 
is by using Bayessian theory. In simple terms this is done 
by forming the so-called Bayes factor[33| Bij, where 


B-ij — 


L(Mj) 

L(Mj) 


( 2 . 12 ) 


and L(Mi) denotes the probability p(D\Mi) (called likeli¬ 
hood for the model Mj) to obtain the data D if the model 
Mi is the true one. Generally, L{Mi) is defined as: 


L{M i )=p{D\M i ) = J da ■ p{a\Mi)Ci{a) (2.13) 


for models with one free parameter and where p(a\Mi) is 
the prior probability for the parameter a. Also, £,(a) is 
the likelihood for the parameter a in the model and 

Ci(a) =e~ x2(a)/2 (2.14) 


In the case that a has flat prior probabilities, that is 
we have no prior information on a besides that it lies in 
some range [a, a + Aa] then p(a\Mi) = and 


i />a+Aa 

L(Mi) = — / da e -x a (o)/ 2 (2.15) 
Aa Ja 


Of course, all this can be generalized for models having 
more than one parameters. 
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The interpretation of the Bayes factor Bij is that [33] 
when 1 < < 3 there is evidence against Mj when 

compared with Mi, but it is only worth a bare mention. 
When 3 < B^ < 20 the evidence against Mj is definite 
but not strong. For 20 < B,j < 150 the evidence is strong 
and for Bij > 150 it is very strong. 

In the next section we will use the Bayes factor B,,j to 
compare various parametrizations to LCDM. 

III. RESULTS-DISCUSSION 

In this section we utilize the p-test to rank 10 rep¬ 
resentative parametrizations of H(z). The comparison 
of the parametrizations considered is shown in Table I. 
This is a representative sample among the many more 
parametrizations we considered but we do not present 
here for clarity reasons. The rank in Table I is accord¬ 
ing to the p-test statistic indicated in the third column. 
A similar rank is obtained using the Xmin/d°f which is 
also sensitive to the number of parameters and shown in 
the last column of Table I. Notice however that the rank 
according to Xmin (shown on the fourth column of Ta¬ 
ble I) which is insensitive to the number of parameters is 
somewhat different. 

The model name initials of the first column correspond 
to the following: ‘OA Var (1)’ is an oscillating ansatz 
H(z) with amplitude decreasing with time like (1 + 2 ) 3 . 
‘OA (2)’ is a similar ansatz but with constant oscillation 
amplitude. ‘LA (3)’ is an ansatz used by Linder [(D[ (but 
proposed earlier in (3f|) and is derived from an equation 
of state parameter of the form 

w{z) = w 0 + Wi—(3.1) 
i + z 

interpolating between two constant values of w(z). ‘P2’ 
is a quadratic polynomial of H{z) also discussed in Ref. 
[36|. ‘Linear’ is an ansatz derived by demanding a linear 
form of the equation of state 

w(z) = wo + w\Z (3-2) 

‘P3’ is a cubic polynomial of H(z). Here we present the 
deeper minimum of % 2 for P3 which shows oscillating 
behavior. We have found another minimum less deep 
which is very similar to the minimum of P2. ‘CA’ is 
an ansatz based on the generalized cardassian cosmology 
01. ‘Quiess’ corresponds to dark energy with constant 
equation of state w. ‘MCG’ is a modified form ;26] of 
the Chaplygin gas J27|. We have tested several variants 
of this ansatz but they all have a fit that is marginally 
distinguishable from LCDM. ‘Brane2’ is an ansatz moti¬ 
vated from brane world cosmology M- It appears in two 
variants which involve a change of sign in the square root 
appearing in H 2 (z) but we find no significant difference 
in the quality of fit between the two variants (‘Branel’ 
and ‘Brane2’); they are both practically indistinguishable 
from LCDM. 


In all cases we have assumed priors corresponding to 
flatness and D 0m = 0.3. For a prior D 0m = 0.25 we have 
found no significant changes in the results of Table I. The 
following are noteworthy features of Table I: 

• The ansatz giving the best fit to the Gold dataset is 
an oscillating parametrization where the amplitude 
of the oscillations decreases like the third power of 
the scale factor. This type of decreasing amplitude 
oscillations can be induced by an oscillating scalar 
held minimally or non-minimally coupled ( e.g. the 
radion jlj, LLil ) ■ With respect to this parametriza¬ 
tion the simplest data consistent parametrization 
(LCDM) is excluded at the 85% confidence level. 
This level of significance is clearly not enough to 
disfavor LCDM but it may provide a hint for inter¬ 
esting better fits to the data. 

• The cubic polynomial parametrization P3 which 
provides a fairly good fit to the Gold dataset (the 
third best Xmin) a ^ so shows indications of H(z) os¬ 
cillations at its best fit despite the fact that no os¬ 
cillating functions are built in the parametrization. 
We have seen a similar behavior in other power law 
parametrizations even though we have not included 
those in Table I to avoid confusion. 

• Physically motivated parametrizations like the 
Chaplygin gas, brane-world models and Cardas¬ 
sian cosmology provide marginally better fits than 
LCDM and are disfavored by our test due to their 
larger number of parameters compared to LCDM. 

• Parametrizations that allow w(z = 0) < —1 and 
crossing of the phantom divide line w = — 1 produce 
significantly better fits to the Gold dataset. 

In Table II, we consider the same parametrizations, as 
in Table I, but now we have marginalized y 2 over D m . 
We have used uniform marginalization in the range D m G 
[0.20, 0.34], We also show the Bayes factor for each model 
against LCDM, calculated with the techniques described 
in the previous section and assuming uniform probability 
distribution of the parameters within their lcr range. 

Also, as expected, after marginalizing x 2 over D m , the 
error bars of the best-fit parameters have increased. For 
example, for the P2 model, in the case of the = 0.3 
prior, we found a\ = —4.16 ± 2.53 and <12 = 1.67 ± 1.03 , 
but after marginalizing a\ = —4.45±2.81 and ci 2 = 1.87± 
1.28. In the same fashion, for the LA model we found in 
the first case wq = — 1.58±0.33 and w\ = 3.29± 1.76 and 
in the latter wq = —1.44 ± 0.44 and w\ = 3.09 ± 1.87. 

In Figs 1 and 2 we show the best fit w(z) and H(z) for 
6 representative parametrizations of Table I. 

Our result for the improved fit of oscillating 
parametrizations is consistent with a previous finding 
by two of us G3 using an earlier Snla dataset 1 2. 'll l24l . 
Clearly, even though LCDM is disfavored at the 85% con¬ 
fidence level with respect to the oscillating parametriza¬ 
tion this is not enough to offer conclusive evidence for the 
existence of oscillations in the cosmological expansion. 
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TABLE I: A comparison of the parametrizations used in the literature. In all cases we have assumed priors 
corresponding to flatness and S2om = 0.3. 


Model 

H 2 (z), (fiom = 0.3) 

p-test 

Xmin 

Best fit parameters 

Xmin/dof 


H 2 (z) = Hi [Oom(l + zf ± 1 - L2om+ 



Oi = 

= 0.13 ±0.07 


OA Var. (1) 

ai(l + 2) 3 [cos(a22 + a37r) — cos(a37r)]] 

0.85 

171.733 

a 2 = 

= 6.83 ± 1.61 

1.115 





03 = 

= 4.57 ±0.07 



H 2 (z) = Hi [flom(l + z) 3 + ai cos(a 2 « + a 3 n)+ 



ai = 

-0.30 ±0.17, 


OA (2) 

(1 — ai cos(a 3 ?r) — flom)] 

0.81 

172.368 

a 2 = 

6.34 ±3.19 

1.119 





a 3 = 

-0.37 ±0.14 



H 2 (z) =Hl[n 0 m(l + z) 3 + 



Wo = 

-1.58 ±0.33 


LA (3) 

±(1 - fi 0 m)(l + z) 3( - 1+w ° +wi) e 3 ' will/(1+z) - 1] ] 

0.79 

173.928 

Wi = 

3.29 ± 1.76 

1.122 


H 2 (z) = H 2 \£lQrn(\ + z)^ + ai(l + z)-\- 



ffll = 

-4.16 ± 2.53 


P2 (4) 

02(1 + z) 2 + (1 — Horn — ai — 02 )} 

0.78 

174.207 

a 2 = 

1.67 ± 1.03 

1.124 


H 2 (z) = Hl[Q 0 m(l + z) 3 + 



Wo = 

-1.40 ±0.25 


Linear (5) 

(1 - fiom)(l + z) 3 ^ 1 + W °- W l) e 3 “l Z ] 

0.75 

174.365 

Wi = 

1.66 ±0.92 

1.125 


H 2 (z) = Hq [Oom(l + ^) 3 + &l(l %) + Cl 2 (l + z) 2 ~h 



«1 = - 

-21.79 ± 17.92 


P3 (6) 

a 3(1 + z) 3 + (1 — Horn — 01 — 02 — 03 )] 

0.74 

173.155 

a 2 = 

14.75 ± 13.07 

1.124 





a 3 = 

-3.13 ±3.07 



H 2 (z) = Hl[Q 0 m(l + z) 3 + 



q = 0 

.0058 ± 0.0014 


CA (7) 

n 0 m(i + z) 3 (( 1 + (nz -1)(1 + z) 3 ^f* - 1 )] 

0.50 

175.758 

n = - 

-119 ±40 

1.134 


H 2 (z) = Hl[Q 0 m(l + z) 3 + 






Quiess (8) 

+ (l-«0m)(l + «) 3(1+U,) ] 

0.15 

177.091 

w = 

-1.02 ±0.11 

1.135 


H 2 (z) = Hl{n 0m (l + z) 3 + 



ai = 

- 0 QQQ2 +0 0008 

- 0 O6O 


GCG (9) 

(1 — f2om)(®l + (1 — Ol)(l + z) 3 ) W ] 

0.031 

177.064 

w = 

18.13 ±4.95 

1.142 


H 2 (z) = Hq [Qom(l + z)3 — -\/oi + 0-2(1 + z) 3 + 



fll = 

29.08 ± 7.30 


Brane 2 (10) 

(1 — ^Om + y/ CLl + 02 )] 

0.027 

177.071 

a 2 = 

-0.097 ±0.459 

1.142 

LCDM 

H 2 (z) = Hl[Q 0 m{l + z) 3 + 1 - Horn] 

- 

177.072 

Horn 

= 0.31 ±0.04 

1.135 
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TABLE II: Like the previous table but now with y 2 marginalized over in the range Q m £ [0.20, 0.34]. Also, 
it contains the Bayes Factor Bij, of each model against LCDM. The Bayes factor is not shown for the Brane 
2 and the P2 parametrizations because the marginalization integrals were ill defined in part of the parameter 
space ( H(z ) becomes complex for some parameter values). 


Model 

H 2 (z) 

Bayes Factor Bij 

Xmin 

Xmin/dof 

OA Var. (1) 

H 2 {z) = Hi [n 0m (l + zf + 1 - + 

ai(l + z) 3 [cos(a2« + a37r) — cos(a37r)]] 

9.03 

176.442 

1.146 

OA (2) 

H 2 (z) = Hi [Dom(l + z) 3 + ai cos(a,2Z + a,3n) + 

(1 — ai cos(a37r) — f2 0 m)] 

2.56 

177.029 

1.150 

LA (3) 

H 2 (z) = Hl{n 0m (i + z) 3 + 

+(1 - n 0 m)(l + z) 3(1+wo+wl) e 3 ' wlll/(1+z) - 1] } 

2.63 

178.676 

1.153 

P2 (4) 

H 2 (z) = Hq {l^Orn (1 + z) 3 + ttl(l + z)~\~ 
a 2 (l + z) 2 + (1 - ^om — ai — a 2 )} 

- 

178.874 

1.154 

Linear (5) 

H 2 (z) = Hl[Qom(l + z) 3 + 

(1 - fiom)(l + z) 3(1+w °- w ^e 3wiZ ] 

2.50 

179.173 

1.156 

P3 (6) 

H 2 (z) = H 2 [Q, 0 m(l + z) 3 + Oi (1 + z) + Cl 2 (1 + z) 2 -\- 
a 3 (l + z) 3 + (1 — ^Om — ai - a 2 — 03 )] 

3.68 

177.859 

1.155 

CA (7) 

H 2 (z) = Hl[fl 0 m(l + z ) 3 + 

«o m (l + Z) 3 (( 1 + (QZ - 1)(1 + z) 3(n - 9) ) * - 1)] 

0.78 

182.308 

1.176 

Quiess (8) 

H 2 (z) = Hl{n 0m (l + z) 3 + 

+(i-n 0m )(i + 2) 3(1+u,) ] 

0.98 

181.614 

1.164 

GCG (9) 

H 2 (z) = Hl[Q 0 m(l + z) 3 + 

(1 — f 2 0 m)(ai + (1 — oi)(l + z) 3 ) w ] 

1.26 

181.394 

1.177 

Brane 2 (10) 

H 2 (z) = Hl[Q 0m (l + z) 3 - v /ai+a 2 (l + «) 3 + 

(1 ^Om ~\~ y/CLl ~\~ CI2)] 

- 

182.326 

1.176 

LCDM 

H 2 (z) = Hl[Q 0 m{ 1 + z) 3 + 1 - Horn] 

1 

182.326 

1.161 
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However, it is recently becoming clear that the pos¬ 
sibility of cosmological expansion oscillations is favored 
by a number of independent additional factors, observa¬ 
tional and theoretical. The theoretical factors include 
the following: 

1. The eternal acceleration predicted by LCDM and 
by most other parametrizations creates problems in 
string theory where asymptotic flatness of space- 
time is required ^]. Theoretical models (using 
eg special quintessence potentials or generalizations 
of the Einstein action and extra dimensions llifl l 
have been proposed (39] to cure this problem by 
cosmological expansion oscillations. Unfortunately, 
these models l39l do not predict crossing of the w = 
— 1 as indicated by the data because they are based 
on minimally coupled scalar fields. 

2. A periodic expansion rate of the universe could re¬ 
solve the cosmic coincidence problem. In an oscil¬ 
lating expansion setup, the present period of ac¬ 
celeration would be nothing surprising; it would 
merely be one of the many other accelerating peri¬ 
ods that also occurred in the past jlij El • 

3. Models with extra dimensions require a mechanism 
to stabilize the size of extra dimensions (the radion 
held). However, since the radion held also cou¬ 
ples to the redshifting matter density (it behaves 
like a Brans-Dicke scalar) its equilibrium point is 
slowly shifting with time and radion oscillations are 
generically excited. The energy density of these os¬ 
cillations is dropping approximately like the third 
power of the scale factor 15]. It is exactly this 
power that gives the best ht of the Snla data! The 
constraints imposed on such low frequency oscilla¬ 
tions by local gravity tests can be relaxed by taking 
into account the increased average matter density 
in the neighborhood of the solar system which can 
signihcantly increase the local effective mass of the 
radion {42 ]. 

The observational factors favoring oscillating expansion 
include the following: 

1. Superimposed oscillations 0] or glitches M in the 
CMB multipole moments Ci can improve the fit 
to the first year Wilkinson Microwave Anisotropies 
Probe (WMAP) data. The corresponding drop in 
the y 2 f° r superposed oscillations was found to be 
Ay 2 ~ 10 which is statistically significant Fhl . 
These oscillations were attributed to new physics 
taking place at the inflationary phase reflecting on 
the produced fluctuations power spectrum. How¬ 
ever, small amplitude oscillations of the cosmologi¬ 
cal expansion could also induce similar effects with 
remnants in the present day expansion rate. 

2. The 128 Mpc periodicity in the galaxy redshift dis¬ 
tribution was first observed in pencil redshift sur¬ 
veys and confirmed more recently by large scale 
redshift surveys O HE ^3, El ■ This periodic dis¬ 
tribution could either be attributed to specific fea¬ 
tures of the matter spectrum (coming eg from the 


acoustic oscillations l49l[ at recombination or from 
exotic physics), or to oscillating cosmological ex¬ 
pansion |50L l51l . Even though the origin of this 
periodicity remains an unresolved issue it should 
be pointed out that a recent study has indicated 
that the probability that such periodicity emerges 
statistically in a LCDM cosmology is less than 0.1% 

WiM 

In conclusion, we have performed an extensive compara¬ 
tive study in parametrization space and have evaluated 
the quality of fit of several H (z) parametrizations to the 
Gold dataset. In comparing the quality of fit we have 
used two statistics: the value of y 2 at its minimum (y min) 
which is independent of the parametrization number of 
parameters and the p-test (P < (Xnnro n )) which evaluates 
the quality of fit with respect to LCDM and depends on 
both xmin an d the number of parametrization parame¬ 
ters n. According to both statistics the best fit to the 
Gold dataset is achieved by an oscillating parametriza¬ 
tion with oscillating amplitude decreasing like the inverse 
cubic power of the scale factor. Even a cubic polynomial 
parametrization is showing oscillating behavior at best 
fit. We stress however that the statistical significance 
of the improved fit of oscillating and other parametriza¬ 
tions we found is at less than 2cr level and is therefore not 
enough to exclude the LCDM parametrization. It does 
provide however a hint towards parametrizations that are 
clearly more probable than LCDM to be realized in na¬ 
ture on the basis of the Gold dataset. 

We have also briefly reviewed the theoretical and ob¬ 
servational evidence that has been accumulating during 
the recent years in support of such oscillations. The par¬ 
ticular type of oscillations that are providing the best 
fit to the equation of state w(z) exhibit crossings of the 
w = — 1 divide line. These crossings can not be repro¬ 
duced in any single field quintessence or phantom model 
and probably require non-minimally coupled theories for 
their physical interpretation. 

The mathematica file with the numeri¬ 
cal analysis of the paper can be found at 
http: / /leandros. physics .uoi.gr/ goldfit .htm 
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